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C/3 . 1 Introduction 

. 

The article is devoted to different aspects of the question: "What can be done with a complex- valued 
matrix by a low rank perturbation?"Q| 

From the works of Thompson [15] we know how the Jordan normal form can be changed by a rank 
pL| ■ k perturbation, see TheoremO Particulary, it follows that one can do everything with a geometrically 

simple spectrum by a rank 1 perturbation, see Corollary [TJ But the situation is quite different if one 
restricts oneself to normal matrices, see Theorem [3] and Corollary [2j We think that Corollary [2] may 
be considered as a finite dimension analogue of the continuous spectrum conservation under compact 
perturbations in Hilbert spaces. For unitary and self-adjoint matrices the inequality of Corollary [5] is 
the only restrictions on "what can be done with a spectrum by a rank k perturbation", see Theorem 2] 
CnJ ■ We don't know if there is an analogue of Theorem 2] for normal matrices. It is worth to mention that 

Corollary [2] for self-adjoint matrices follows from Cauchy interlacing theorem [2]. Theorem[4jis related 
with the converse Cauchy interlacing theorem [6] . 

The spectrum of Hi + H2 with known spectra of self-adjoint matrices Hi and H2 is studied a 
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(y-) ■ lot, see [9] and the bibliography therein. Although the complete set of restrictions on the spectrum 



H\ +H2 known in this situation, we are not sure that there is an easy proof of Theorem fusing results 
of 0. 

Although Theorem [2] should be known (see, for example, [13], where Theorem [2] formulated in one 
direction), we will give a proof here, manly because our proof falls in a general framework , which is 
also used in the proof of TheoremHJ Let us describe the framework. The set C„ xn of all complex nxn- 
matrices (set of self-adjoint matrices) we equip with the arithmetic distance, d(A, B) = rank(A — B) 
(see [3]). The arithmetic distance is geodesic for these cases. The spectral properties of matrices, 
such as Weyr characteristics and spectra (multiset) also may be considered as a metric spaces with 
distance, related to the arithmetic distance on matrices, see Section [2] These distances also turn out 
to be geodesic. Then we prove Theorem [2] (Theorem |4j for rank(A — B) = 1 and the general results 
will follow from Proposition [TJ 

Proposition 1. Let X and Y be geodesic metric spaces, let (x) denote the closed ball of radius n 
around x in X , let (f> : X -> Y be such that <f>{0?{x)) = Of((p(x)) for all x G X. Then 4>(0* (x)) = 
0^(4>(x)) for any n G N and x £ X . 

Proof. The proof is by induction. For n — 1 there is nothing to prove. Step n — ► n + 1: It follows 
that 0% +1 (x) = [J Of (z) {X is geodesic), then 

z£0%(x) 

d>(0^ +1 (x)) = |J 0(Of(*))= (J °i>)= U OY(z) = Ol +1 (d>(x)). 

zeO^(x) z£<f>(OX(x)) zdOl(4>(x)) 
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work was partially supported by CONACyT grant SEP-25750, and PROMEP grant UASLP-CA-21. 
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In Section [5] and Section[6]we use the normalized arithmetic distance d r (A, B) = rank ( A ~ B ) ; w here 
n is the size of the matrices. We are interested in the following questions: " Suppose that matrices 
almost satisfy some equations (in the sense of d r (-, ■)). If close to that matrices there exist matrices 
satisfying the equations (uniformly with respect to n)?" We manage to answer only the following: 
close to an almost unitary (self-adjoint) matrix there exists an unitary (self-adjoint) matrix. We do 
not know if the same is true for normal matrices. (This question has the affirmative answer for norm 
distance dn(A, B) = \\A — B\\, see [10]. It is equivalent to the following: "close to any pair of almost 
commuting self-adjoint matrices there exists a pair of commuting self-adjoint matrices (with respect 
to the distance dn(-,-)). It is interesting that there are almost commuting (with respect to dn(-,-)) 
matrices, close to which there are no commuting matrices, 4, 5[[16]). The similar question have been 
studied for operators in Hilbert spaces (Calkin algebras, |7]). In Hubert spaces the operator a is 
called to be essentially normal iff aa* — a*a is a compact operator. In contrast with Theorem[7J there 
exists an essentially unitary operator which is not a compact perturbation of an unitary operator 
(just infinite 0- Jordan cell). There is a complete characterization of compact perturbations of normal 
operators, see [7] and the bibliography therein. Let us return to almost commuting matrices with 
respect normalized arithmetic distance d r . In Section [6] we show that for any A 6 C nxn with simple 
spectrum there exists an almost commuting matrix, which is far from each commuting with A matrix. 
The similar problem for the pairs of almost commuting matrices, as far as we know, is open. Precisely, 
if for any e > there exists 5 > such that for any A, B, d r (AB, BA) < e there exists (A, B) satisfying 
AB = BA and d r (A, A),d r (B,B) < 5 (5 does not depends on size of the matrices). 

We think that low rank perturbations of matrices may related to sofic groups. The following 
question seems to be interesting from this point of view (although it seems to goes beyond the scope 
of the present article). One can show that all solutions of equation C~ 1 A~ 1 CAC~ 1 AC = A 2 in finite 
unitary matrices are trivial in A (A = E). On the other hand, it is true that for any e > there exist 
A,C, d(A, E) = d{C,E) = 1 and d(C~ 1 A~ 1 CAC~ 1 AC, A 2 ) < e. If the above assertion is true with 
additional requirements C 4 = 1? If not, it gives an example of non-sofic group. 

Note. All linear spaces are supposed to be finite dimensional in the rest of the article. C„ xn will 
denote the set of all complex n x n-matrices, N = {0,1,2,...}. 

2 Some discrete geodesic spaces. 
2.1 Arithmetic distance on C nxn 

Lemma 1. The arithmetics distance rank(A — B) is geodesic on 

• Set of all n x n matrices. 

• Set of all self-adjoint n x n matrices. 

• Set of all unitary n x n matrices. 

Proof. It is clear that a rank k matrix (self-adjoint matrix) may be represented as sum of k matri- 
ces (self-adjoint matrices) of rank 1. The first two items follow from the fact that set of matrices 
(self-adjoint matrices) is closed with respect to summation. For unitary matrices. Let rank{U\ — 
U2) = fc, or, the same, rank(E — fTf U2) = k. It means that, in a proper basis, U{~ U% — 
diag(Xi, A2, Afe, 1, 1, 1). Now the sequence U\, Urdiag(\\, 1, 1, 1), Ui-diag(Xi, A2, 1, 1, 1)...U%- 
diag(Xi 7 X 2 , Afe, 1, 1, 1) = U2 give us the geodesic needed. □ 

Remark 1. The methods used in the above proof are not applied for normal matrices - the set of 
normal matrices is not closed with respect neither summation nor multiplication. In fact, an example 
from [5] hints that arithmetic distance might be non geodesic on the set of normal matrices. 
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Proposition 2. Let <f)(x) = (ax + b) (cx + d) be a Mobius transformation of C nxn , defined on A, 
B. Then rank(A - B) = rank((/)(A) - <j)(B)). 

Proof. A Mobius transformation is a composition of linear transformations A — ► aA + b (a, & € C) 
and taking inverse A — > A^ 1 . Those transformations (if defined) clearly conserve arithmetic distance, 
for example, rank(A~ x - B^ 1 ) = rank(A^ 1 (B - A)B~ r ) = rank(A - B). (A^ 1 and B^ 1 is of full 
rank.) □ 



2.2 Distance on the spaces of the Weyr characteristics. 

Having in mind the Weyr characteristics of complex matrices (see below) , we introduce the spaces Ss„ 
of the Weyr characteristics. Where 9„ is the space of functions Z+ x C — > N, (i,X) — > %(A) such 
that 

• ?7i(A) 7^ for finitely many (i, A) only, and ^ ^2 T)i(\) = n. 

AeCieN 

• ViW > Vi+iW- 

On 3 ra define a metric d(n,fi) — max{|?7i(A) — fii(X)\}. First of all let us note that d(-,-) is indeed 

(i,A) 

a metric. Trivially, d(rj,fx) = implies rj = /i and d(-, •) satisfies triangle inequality as supremum 
(maximum) of semimetrics. It is clear, that d(fi, v) is also well defined for fi and v in different spaces 
of Weyr characteristics (for different n) . We will need the following 

Proposition 3. Let fi G 5 m and n > m. Then there exists v G 3„ such that for any rj G 5„, the 
inequality d(fi,rj) > d(y,rf) holds. 

Proof. We can do as follows. Let /x,(Ao) ^ and /Xj + i(Ao) = 0. We can take i/£ + i(Ao) = z^+2(Ao) = 
... = !/j +n _ m (Ao) = 1 and /ij(A) = fj(A) for all other pairs (J, A). □ 

Proposition 4. Q'n are geodesic metric spaces. 

Proof. Let 77,// G 3„ and d(r],(J,) = k > 1 it is enough to find 2/ E 3 n such that either d(rj,v) = 1 
and d(v, n) = k — 1, or d(ri,v) = k — 1 and d(v, fi) = 1, moreover, by Proposition [3] it is enough 
to find t/ G 3 m for m < n. Let 5+ = {(j, A) e Z+ x C | ^(A) - ^(A) = fc} and S_ = {(j, A) G 
Z + x C I ^j(A) — Hj(X) — ~k). Suppose, that \S+\ > \S-\ (if not, we can change rj «-» fi). Now let 

f 7fc(A)-l if (*,A)eS+ 
i/i(A) = < ??i(A) + l if (i,A) € S_ 

[ jji(A) if (i,A)^ + U5_ 

We have to show that v G Ss m for m = n — \S+ \ + \S-\. It is enough to show that ja, +1 (A) < fj(A). 
Suppose contrary, that tA,- + i(A) > Vj(\). There are three possibility: 

a) %+i(A) = r]j(X), (j, A) G S+ and (j + 1, A) g" 5+, but then fc > %+i(A) - Mj+i(A) > %+i(A) - 

Mi (A) = *?i(A) — Mj(A) = fc, a contradiction. 

b) %+i(A) = »y(A), + A) G 5_ and (j, A) £ S_, but then -k < (A) - ^ (A) < ^(A) - Mj+i(A) = 

^3+1 (A) — /ij+i(A) = — fc, a contradiction. 

c ) ? /j+i(A) = VjW - 1, (j)A) G 5+ and (j + 1,A) G 5_, but then -k = Vj+iW - Mj+i(A) > 

r; J+ i(A) — i-J.j(X) = r)j(X) — (J,j(\) — 1 = k — 1, so — fc > A; — 1 and 1/2 > k, a contradiction with 
k > 1. 

Now, by construction, d(j/, 77) = 1 and d(z/, /1) = fe — 1. □ 
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2.3 Distances dc and dc on finite multisets of the complex numbers. 

Multisets and operations. The language of multisets is very convenient to deal with spectrums. We 
will need only finite multisets. For a multiset A let set(A) denote the set of elements of A, forgetting 
multiplicity. It is clear that a muiltiset maybe considered as the multiplicity function xa '■ set(A) — > N, 
for any x A we will suppose Xa( x ) = 0. (For all cases, considered here, set (A) C C, so we can 
consider xa : C — > N = {0, 1, 2...}.) As far as the authors aware, there are several generalizations of 
the set-theoretical operations to multisets. We will need the following operations: 

• Difference of two multiset A \ B, Xa\b{ x ) = niax{0, x(-4) — x(B))}- 

• Intersection A fl X of a set X and a multiset A, 



XAnx(x) 



Union A&B, xa<±>b(x) = xa(x) + Xb(x). 



Xa{x), if x G X 
0, if i^I 



Let S(a,r) = {x G C : \x — a\ < r} and S = {S(a,r) \ a £ C, re R+} be the set of circles. For 
A,B Cm C let dc(A,B) — max{||.4n 5| — ISnS 1 !!}. Let as extend S to S which include interior of 

complements of circles and semiplains: S = S U {{x G C : \x — a\ > r} | a G C, r <E R + } U {{x G 

C : Im(^) > 0} | a, b G C}; and introduce new metric dc: dc(A, B) = max{||^ D S\ - \B D S\\}. 

ses 

Proposition 5. • dc and dc are metrics on the finite multisubsets of C 

• dc{A,B) = dc{A\B,B\A), dc{A, B) = dc(A\B,B\A). 

• // \A\ = \B\, then dc(A,B) = dc(A,B). 

Proof. • The same as for the spaces of Weyr characteristics. 

. LetEs(AB) = \AHS\-\BnS\ = E (xa(x)-Xb(x)). Then J2 S (A\B, B\A) = E (max{0, X a(x)~ 

x£S x<£S 

Xb(x)} - max{0, xb(x) - Xa(x)}) = E (Xa(x) - Xb{x)). Now the item follows by definition of 

xes 

dc (dc). 

• First of all, due to A and B arc finite multisets, for any semiplain p we can find a circle c such 
that E P (A ^) = S c ("4' ^O- Al so f° r an y closed circle c c there exists an open circle c such that 
E Cc (-4j<B) = E Co (-^'^)- Now, under assumption of the item ^2 S {A 7 B) = — Ec\s(^'^) an< ^ 
the result follows. 

□ 

Proposition 6. Let 4>(x) = be a Mobius transformation ofC, defined on set(A)L) set(B). Then 
dc(A, B) = dc{4>(A),4>{B)). 

Proof. A Mobius transformation defines a bijection of S. □ 

We don't know if the metric dc is geodesic on the multisets with fixed cardinality, but its restriction 
on any circle o line is: 

Proposition 7. Let I C C be a circumference or a straight line. Let A,B Cm h \A\ = \B\ = n and 
dc(A,B) = k > 2. Then there exists C Cm I, \C\ = n such that dc(A,C) = 1 and dc(C,B) = k — 1. 
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Proof. By Proposition [6j it is enough to proof it for the unit circle. Let us start with the case when 
set(A) n set(B) = 0. Let T = set(A) U set(B) C C 1 . Let |L| = r. We will cyclically anticlockwise 
order P = {70, 71, J r -i} by elements of Z r . To construct C we move each element of A to the next 
element in L, precisely, set(C) C L and 

Xcili) = max{0,x^(7 4 ) - 1} + X S et(A)(li-l)> 

the other words 

!Xa{h) - 1 if Ji 6 sef(.A) and 7^ ^ set(^) 
1 if 7^ ^ sef(.A) and 7^1 6 se£(U) 

X-4(7«) f° r t ne other cases 

We will check that C satisfies our needs. For x, y € T let [x, y] denote the closed segment of C 1 , starting 
from x and going anticlockwise to y (so [a;,y] U [y,x] — C 1 ). It is clear that for X, Y Cm L one has 
dc(X,Y) = max{||Xn [a,/3]\ - \Y n [a,0}\\ : a, (3 e set(X) U sei(F)}. Denote by E[ a /3]( X > y ) = 

|xn[a,^]J-|yn[a,/3]|. 

Now, dc(A,C) = dc(A \ C,C\ A) = 1, for A \ C and C \ A are interlacing sets on C 1 . Suppose 
further, dc{C, B) — dc(C \B 7 B\C) > k, then there exists [7*, 7j] such that either 

1- E[ 74l7j] (C\B,B\C) = dc(C \B,B\C)>k 
or 

2- E [74)7j] (C\B,B\C)<-fe. 

In the first case we may assume that 7^, 6 C \ B and 74-1, 7j+i ^ C \ B. Now, changing interval if 
necessary (i n — i — 1 and (or) j™ = j — 1) we may, keeping E[ 7i 7 -] (C\B,B\C), achieve that 74, 7,- G „4 
and 7i_i,7j +1 A Then Er 74t7J ](-4iB) = E[ 7l , 7j] (C \ B,B \C) + 1 > k + 1, a contradiction. The 
second case may be considered similarly. 

If set(A) n set(B) ^ then we can find C for .4 \ B and B\A and then take C = C fell, where 
X = .A\(.A\£)=£\(S\.4). □ 



3 On spectrum of low rank perturbations. 

Theorem 1 (Thompson). Let n x n matrix A over a field F have similarity invariants h n (A) \ 
h n —i(A) I ... I hi(A). Then: as column n tuple x and row n-tuple y range over all vectors entries in 
F , the similarity invariants assumed by the matrix 

B = A + xy 

are precisely the monic polynomials h n {B) . . . h\{B) over F for which degree{h\(B) ■ ■ ■ h„ (B)) = n 
and 

h n {B) I h n -x(A) I hn-2(B) I h n - 3 (A) I 
K{A) I K-i{B) I K- 2 {A) I K-z{B) I .... 

We are going to reformulate Theorem Q] for the field C using Weyr characteristic. 

Let r] m (A, A) denote the number of A- Jordan blocks in A of size greater or equal to m (m G N). 

r} m (A, A) = dim Ker(XE - A) m - dim Ker{\E - A)" 1 ' 1 

This sequence of numbers r)i(A, A), . . . f] q (A, A) is called the Weyr characteristic for the eigenvalue A 
of matrix A, see [13]. 

Theorem 2. Let A e C„ x „ with Weyr invariants rj m (A, A). Then as R ranges over all nxn complex 
matrices of rank less o equal k, the Weyr invariants assumed by the matrix B = A + R are precisely 
those, that satisfy both of the following conditions: 
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q n -i(B,\) 
Qn(B,X) 



Figure 1: Relation between Ferrers diagrams of (B, A) and (A, A) 



For any A € C and any m£N 



| T] m (A,X)-r] m (B,X) |< fc. 



• L Vm{B, A) = n. 

AeCmiEN 

Proof. It is enough to prove the theorem for k = 1. Indeed, assume that the theorem is valid for 
k = 1. Then we may consider the Weyr characteristics as a map ip : C nxn —* 3„, which satisfies 
Proposition [TJ So, Theorem [2] follows, though it states that tp(Ok(A)) — Ok(*p(A)). 

Let us prove Theorem [2] for k = 1. For a given eigenvalue A € sp(A) the sequence of numbers 

qi(A,X)>q 2 {A,X) > 

corresponding to the sizes of the A- Jordan blocks in the Jordan normal form of A are know as the 
Segre characteristics of A relative to A, [14] . 

The similarity invariant factors of A G C nx „ are sequence of monic polynomials in x, h n (A) \ 
/i„_i(A) | h n - 2 (A)... | It is known that h,(A) = ]J X (X - x ) q '^ x) , where A G sp(A) and 

qi(A, A) is a Segre characteristic corresponding to A. So, if rank(A — B) = 1, by Thompson's theorem 
[T]one has 

91 OB, A) >q 2 {A,X) >q 3 (B,X) >..., 
qi{AX) >q 2 {B,X)>q 3 {A,X) > .... 

As Weyr characteristic is the conjugate partition of Segre characteristic, we can use Ferrers diagram 
to compute w m {B, A) (see, [T3]) as the number of points of column m in Ferrer diagram of the Segre 
characteristic of B relatively to A (by short the Ferrer diagram of (B, A)). Precisely, the Ferrer diagram 
for q(B, A) is the set Fg = G Z + x Z + | j < qi(B, A)}, see Figured] (the numbering from top 

to bottom and left to right). The Weyr characteristics is related with Ferrers diagram by the formula 
r]j(B,X) = \{(x,y) G Fg \ x — j}\. From inequalities Jl]) we have that qi(B.X) > qi + i(A, A). This 
inequality is equivalent to the statement Vi ^ 1 (i, j) G F% ~ > (i — G Fg (Figure [T]), which is 
equivalent to the fact that rjj(B,X) > i]j(A 1 X) — 1. In similar form, from inequalities ([T]) we can 
observe that qi(A, A) > qi + i(B, A) with 1 < i < n — 1. Therefore, r/j(A, A) > r]j(B,X) — 1, and the 
theorem follows. 

□ 

Corollary 1. If the geometric multiplicity of any eigenvalue X of A (number of X- Jordan cells)is 1, 
then for any multiset M of size n there is a rank 1 matrix B such that sp{A + B) = M . 



(1) 



4 Case of normal matrices. 

We will say that the vector x is an a-eigenvector if Ax — ax. We will denote by 1Z(A, A, e) the 
space generated by all the a-eigenvectors of A with |A — a\ < e. For the case of normal matrices, the 



G 



following theorem shows that the difference between the dimension of TZ(A, A, e) and the dimension 
1Z(B, A, e) is bounded by the rank of the difference matrix A — B. 

Theorem 3. If A and B are normal matrices, then for any A, and for any e > 0, 
| dim(n(A, A, e)) - dim(K(B, A, e)) |< rank(A - B) 
Let X 1 - be the orthogonal complement of subspace X and Px be an orthogonal projection on X. 

Lemma 2. Let N : L — » L be a normal operator, and X be a subspace of L such that \\(N —X)x\\ < e||x|| 
for any x G X , then (we will write e) for 1Z(N, A, e) ) 

1- Ptz(\^)X ± for any x G X, x ^ 0. 

2- \\Piz(x,ae)x\\ > y/l- J||a;|| for any xeX. 
3. dim(R(\,e)) > dim(X) 

Proof. It is clear that (1) implies (3). 

(1) Let ei , e 2 , e n be a diagonal orthonormal basis for N, and Ai, A2, ...A n corresponding eigen- 

n 

values (Nei = A^). Let x = o^ei + ot2e 2 + ... + a n e n G X and ||x|| = \ a i\ 2 = 1- Let x = x\ + x 2 

i=l 

n 

where xi G K(\,e) and x 2 G Tl^{\,e). Now, \\{N - \)x\\ 2 = l^| 2 |Ai - A| 2 < e 2 implies that 

i=i 

J2 |a:.i I 2 < 1. So, Pn{\,e){x) =xi= a i e i 0- 

t I |Aj-A|>e i I |Ai-A|<e 

(2) Similarly, 

i«,i 2 <4 and E m 2 >(i-4)< 

i I I Ai — A|>ae i | | A* — A|<ae 



BO, ||P w(A , ae) x||>^l-^||x|| □ 

Now we a ready to prove Theorem [3] Let rank(A — B) = r then there exists X = TZ(A,X,e) D 
ker(A-B) with dim(X) > dim(1Z(A, A, e)) - r and A\ x = B\ x , so ||(S-A)|x|| < e and, by LemmaEl 
dim(lZ(B, A, e)) > dim(lZ(A, A, e)) — r. By symmetry, we get Theorem[3l 

Theorem [3] implies that any circle in complex plain, containing m spectral points of A (B) should 
contain at least m — k spectral points of B (A). So, we have 

Corollary 2. If A and B are normal matrices then dc(sp(A) , sp(B)) < rank(A — B). 

Proof. It is just a reformulation of Theorem [3) □ 

If the condition of Corollary [2] describes all accessible by rank k perturbation spectra? We are 
going to show that the answer is "yes" for self-adjoint and unitary matrices. 

Theorem 4. Let A be a self-adjoint (unitary) n x n-matrix. Let B Cm R- (B Cm G 1 ), \B\ = n. Then 
there exists self-adjoint (unitary) matrix B such that sp(B) — B and rank(A — B) = dc(sp{A),B). 

In fact the following, more general result is valid: 

Theorem 5. Let I C C be a circumference or straight line. Let A be a normal n x n-matrix, 
sp(A) Cm I- Let B (1m I, \B\ = n. Then exists a normal matrix B such that sp(B) = B and 
rank(A - B) = dc(sp(A),B). 
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• It is enough to prove Theorem [5] for self adjoin matrices. Indeed, let sp(A), B Cm I, \B\ — n for 
a circle (line) I C C. Then there exists a Mobius transformation <f>, defined on sp(A) UB, which 
map I to the real line. Then <p(A) is a self-adjoint matrix and we can apply Theorem 0] to (f>(A) 
and <t>{B) to find B with sp{B) = <j>{B) and rank{cj>(A) — B) = dc{<f>{sp{A)) , 4>{B)) . Now take 
B = </> _1 (i?) and results follows, for the Mobius transformations conserve arithmetic distance 
on C nxn and the distance dc on multisets (Proposition [2] and Proposition [6]). 

• It is enough to prove Theorem [5] for dc(sp(A), B) — 1 and the rest will follow from Proposition [TJ 
Proposition [7] and Lemma [TJ 

• Also w.l.g. we may assume that set(sp(A)) n set[B) = 0. For if X = sp(A) \ (sp(A) \ B) we 
can write A = A\ © A 2 with sp{A\) — X and sp(A 2 ) — sp(A) \ X. We can find B 2 with 
sp(B 2 ) =B\X and rank(A 2 - B 2 ) = 1. Now, take B = A x ® B 2 . 

• Let A,B Cm R, set(A) n set(B) = 0, = \B\ and dc(A,B) = 1. Then, in fact, A and B are 
interlacing sets. It means that if A = {ai,ct2, ...,a n }, B = {Pi,P 2 , then a>i < f3\ < a 2 < 
02 < ■■■ or Pi < a.\ < (3 2 < a 2 < .... 

So, we need to prove only 

Lemma 3. Let A 6 C nx „ be a self-adjoint matrix with a simple spectrum. Let 6cR with \B\ — 
n. Lf sp(A) and B are interlacing then there exists a self-adjoint matrix B with sp(B) = B and 
rank(A - B) = 1. 

Let sp(A) = A = {oti, a 2 , a n } and B = {[3i, j3%, /3 n }. As A and B can be put in diagonal 
normal form A — diag{ai,a 2 , ...,a n ) and B — diag(/3i, (3 2 , (3 n ) by unitary transformations and 
unitary transformations map (by conjugation) self-adjoint matrices to self adjoint matrices, Lemma [3] 
is equivalent to the fact that under our assumptions on A and B the equation 

AX — XB = R (2) 

has a solution in (X, R) for unitary X and R of rank 1. Before solving Eq[2]let us introduce some 
notations and prove a proposition. For a finite A C R let -Pa(A) = Y\ a£ ^{^~ a )- Let A and B be finite 
subsets of R of equal cardinality. It follows from interpolation that there exist unique x : A — ► A, 
such that 

P, b = Pa - / ] x a P A \ {a} (3) 

(we write x a not a;(a)). Moreover, x a = Ps{a)/Pj l \s a \{ot). Studying signs of Pb(o.) and P^\[ a y(a) 
we trivially get 

Proposition 8. If A and B interlacing, then all x a in Eq^have the same sign. 

Remark 2. In fact inverse of this proposition is also valid, see Lemma 1.20 of [8]. 

Now let us go back to solutions of EqO For R = {r^} fixed the equation has the unique solution 
X = {xij}, with Xij — - . Now, suppose that rank(R) — 1 or, the same r%j = ytzj for some 
y, z £ C n , y, z =/= 0. When the matrix X is unitary? When its columns (rows) are orthonormal, or 

It follows that Zj ^ for all j = l,...,n, changing rows by columns we get the same for y. So, the 
difficult part is to guarantee that l.h.s. of Eq. [4] is for j ^ k. Putting equality 

\y*\ 2 i ( kl 2 y 



(Oti - Pj){0Li - Pk) Pj - Pk OLi - pj Oii-Pk 
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into Eq[4]and multiplying it by PA{{3k)PA((3j)(flj — h)/{zjz1) we get (after some elementary trans- 
formations): 

n n 
z=l i=l 

which imply that there exists c, such that 

n 

Y J \yi\ 2p A\{o H }{Pj) = cPA{[3j), 

1=1 

for all j. Or the same, B is the set of roots of polynomial 

n 

$(A) = J] Iwl^Uo*}^-) - cP4(/3i), 
»=i 

so \yi\ 2 = Pg(aj)/ 'cP4\{ Qi }(a,). Choosing c = 1 or c = —1, we get, by Proposition [51 that \yi\ 2 is 
well defined. Now, take, for example, g/j = From Eq0]for j = k we can find \zj\ 2 . Then, taking 
Zj = \zj\ we get needed solution of Eq[2j 

5 Almost unitary operators are near unitary operators with 
respect to normalized arithmetic distance 

For A, B e C„ X n, let d r (A, B) be the normalized arithmetic distance: 

d r (A,B)= rank{A - B) 
n 

The matrix A is called an a-self-adjoint matrix if d r (A, A*) — a, where A* denotes the adjoint of 
A. The matrix A is called an a- unitary matrix if d r (A* A, E) = a 

The following theorems says that "near" to any a-self-adjoint matrix there exists a self-adjoint 
matrix S, and that "near" to any a- unitary matrix there exists an unitary matrix U (for small a). 

Theorem 6. For any A S C nxn there exists a self-adjoint matrix S (S = S* ) such that d r (A,S) < 
d r (A,A*). 

Proof. Take S = | (A + A* ). □ 

Theorem 7. For any A £ C nxn there exists a unitary matrix U (U*U = E), such that d r (A, U) < 
d r (A*A, E). 



The good illustrations for this theorem are 0- Jordan cells: 
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but the matrix 

/0 ••• ••• f\ 
f ••• ••• 



••• ••• 
■•• 10 ••• 



\0 ■•• •■■ 0/ 
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is unitary. 



Proof. Let rank (A* A — E) = r, so there exist subspace Icl, dim(X) — n — r such that ^4*A|x = 
E\x- Consider A\x ■ X — > Y = A(X). Under assumptions of the theorem A*{Y) = X, it follows 
that = A*\y : Y — > X, so A\x ■ X — > Y is an unitary operator. Choose any unitary operator 

B : X 1 - -> Y 1 - (B*B = E x ±). Then U = A\ x © B proves the theorem. □ 

It is not clear if this proof could be adapted for normal matrices - unitary operator from an unitary 
space to another unitary space is well defined, but how to define normal operators between different 
unitary spaces...? 

Question: If we define a-normal matrices in similar form to self-adjoint and unitary matrices, the 
equivalent of theorems [5] and [7] are true for normal matrices? 



6 Almost commuting matrices 

Theorem 8. For every 4 < n G N and every A £ C nX n with simple spectrum there exists X £ C n xn 
such that d r (AX,XA) < 2/n and for any matrix B, commuting with A, d r {B,X) > h. 

Before starting the proof of the theorem we need some facts. 

Proposition 9. Let {Ai, A&} and {ai, ...,ak} be two disjoint sets, then the matrix M = [xij] with 



x 



is nonsingular 



% 3 — ai-Xj 

Proof. The matrix M has the form 

M = 



/ 1 

qi t Ai 

V 



ai — A2 
1 1 

Oi2—\l Ot2 — A2 

1 1 

Oik— Al 



a-2—\k 
1 



otk— A2 



oik — 



/ 



Let P(X) = (A — Ai)(A — A2) • • ■ (A — Afe). Multiply each row j of the matrix M by P(otj) we obtain a 
matrix of the form 

( Pi{ai) P 2 {a 1 ) ••• P fe (ai) \ 

M = 



Pi(a 2 ) P 2 (a 2 ) 
V Pi{a k ) P 2 (a k ) 



Pk(a 2 ) 
Pk{a k ) J 



where P 3 (\) = ^ 



This matrix will be nonsingular if and only if matrix M is nonsingular. We 

will prove that matrix M is nonsingular showing that the following system of linear equations has a 
unique solution: 

( P 1 (a 1 ) P 2 (ai) ••• P k (ai) \ ( ai \ ( bl \ 
Pi(a 2 ) P 2 (a 2 ) ■•■ P k (a 2 ) «2 



\ Pi(a k ) P 2 {a k ) ■■■ P k (a k ) J \ ak J \ b k J 

So we have to solve the system JD i=1 a>iPi(&j) — bj. Consider the polynomial ^(A) = Yli=i a iPiWt 
note that ^(Ai) = a^P^Ai) because Pi(Xj) = for i ^ j. We have k points bj, therefore, we can use 
Lagrange interpolation to find the unique polynomial ^(A) of degree k — 1 such that &{ctj) — bj, and 
then we can compute the values a, = -jrp^- (Pj(Aj) ^ for all Ai are different). □ 



Now we are ready to prove Theorem |8l 
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Proof. Consider the matrix equation 



AX — XA = {cij}, 



(5) 



with Cij — i + jmod2. This matrix has the following form 



AX — XA 



1 1 

1 

1 1 
1 



1 + n mod2 \ 

n mod2 
1 + n mod2 

n mod2 



V n + 1 mod2 







Let A be a diagonal matrix with simple spectrum A — diag(X%, \%, X n ). One matrix X that 
satisfies Eq. [5] is X = {xij} with 

for i =/= j 



Xi — Xj 







for i = j 



Every matrix B that commute with A should be necessarily a diagonal matrix B — diag(bi, 62, b ri 
then X — B = {x* 3 } with 

I x C '\ for i 7^ j 
-h for i = j 



1 



If we delete from X — B the odd columns and the even rows, we obtain a submatrix X' of size 
L§J x Lf J o f the form 

IV 1 = 

A* - A* 

with i* — 2i — 1, j* = 2j. By proposition [9] this matrix is nonsingular. Therefore, rank(X — B) > ^, 
and we obtain d r (X, B) > h. □ 
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